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ABSTRACT 

We investigate the prospects for characterizing extrasolar giant planets by measuring planetary oblateness 
from transit photometry and inferring planetary rotational periods. The rotation rates of planets in the solar 
system vary widely, reflecting the planets' diverse formational and evolutionary histories. A measured oblateness, 
assumed composition, and equation of state yields a rotation rate from the Darwin-Radau relation. The lightcurve 
of a transiting oblate planet should differ significantly from that of a spherical one with the same cross-sectional 
area under identical stellar and orbital conditions. However, if the stellar and orbital parameters are not known a 
priori, fitting for them allows changes in the stellar radius, planetary radius, impact parameter, and stellar limb 
darkening parameters to mimic the transit signature of an oblate planet, diminishing the oblateness signature. 
Thus even if HD209458b had an oblateness of 0.1 instead of our predicted 0.003, it would introduce a detectable 
departure from a model spherical lightcurve at the level of only one part in 10^. Planets with nonzero obliquity 
break this degeneracy because their ingress lightcurve is asymmetric relative to that from egress, and their best- 
case detectability is of order 10^**. However, the measured rotation rate for these objects is non-unique due to 
degeneracy between obliquity and oblateness and the unknown component of obliquity along the line of sight. 
Detectability of oblateness is maximized for planets transiting near an impact parameter of 0.7 regardless of 
obliquity. Future measurements of oblateness will be challenging because the signal is near the photometric limits 
of current hardware and inherent stellar noise levels. 

Subject headings: occultations — planets and satellites: general — planets and satellites: individual (HD209458b) 
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1. INTRODUCTION 

Discovery of a transiting extrasolar planet, HD209458b 
(Charbonneau et al. 2000; Henry et al. 2000), has pro- 
.vided one mechanism for researchers to move beyond 
'the discovery and into the characterization of extraso- 
lar planets. Precise Hubble Space Telescope measure- 
ments of HD209458b's transit lightcurve revealed the ra- 
dius (1.347 ± O.OeORjup) and orbital inclination (86.68° ± 
0.14°, from impact parameter 0.508) of the planet, which, 
along with radial velocity measurements, unambiguously 
determine the planet's mass (0.69±0.05Mjup) and density 
(0.35g cm-3; Brown et al. 2001). Of the over 100 ex- 
trasolar planets discovered so far, this large radius makes 
HD209458b the only one empirically determined to be a 
gas giant. 

Knowledge of a planet's cross-sectional area provides a 
zeroeth order determination of its geometry, and the HST 
photometry is precise enough to constrain the shape of 
HD209458b to be rounded to first order. While a planet 
transits the limb of its star, the rate of decrease in appar- 
ent stellar brightness is related to the rate of increase in 
stellar surface area covered by the planet in the same time 
interval. We investigate whether it is possible to use this 
information to determine the shape of the planet to second 
order. 



Rotation causes a planet's shape to be flattened, or 
oblate, by reducing the effective gravitational acceleration 
at the equator (as a result of centrifugal acceleration) and 
by redistributing mass within the planet (which changes 
the gravity field). Oblateness, /, is defined as a function 
of the equatorial radius (Req) and the polar radius (i?p): 
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For Jupiter and Saturn, high rotation rates and low densi- 
ties result in highly oblate planets: fjupiter = 0.06487 and 
f Saturn = 0.09796, Compared to f Earth = 0.00335. 

An earlier investigation of the detectability of oblateness 
in transiting extrasolar planets was published by Seager & 
Hui (2002). Our work represents an improvement over 
Seager & Hui (2002) in the use of model fits to compare 
oblate and spherical planet transits, the use of the Darwin- 
Radau relation to associate oblateness and rotation, and 
a thorough investigation of the degeneracies involved in 
fitting a transit lightcurve. 

In this paper, we investigate the reasons for measuring 
planetary rotation rates, the relationship between rotation 
rate and oblateness for extrasolar giant planets, the effect 
of oblateness on transit lightcurves, and the prospects for 
determinating the oblateness of a planet from transit pho- 
tometry. 
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Fig. 1. — Spin to orbital mean motion ratios for tidally 

evolved fluid planets in equilibrium. The solid line 
represents the ratio as calculated under the assump- 
tion of a frequency- indenpendant tidal dissipation fac- 
tor Q, and the dashed line is calculated assuming Q oc 
frequency^^. Under these assumptions, extremely eccen- 
tric planet HD80606b (e = 0.93) would, if allowed to come 
to tidal equilibrium, rotate over 90 times faster than its 
mean orbital motion! 
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Fig. 2. — The moment of inertia coefficient, C, as a func- 
tion of planet mass for hypothetical generic 1.0 Rjup ex- 
trasolar giant planets. Extrasolar giant planets may or 
may not possess rocky cores depending on their formation 
mechanism, so we plot C for both no core (upper curve) 
and an assumed 20 Mq core (lower curve). 



2. EXOPLANETS AND ROTATION 

While knowledge of a planet's mass and radius provides 
information regarding composition and thermal evolution, 
measurements of rotation and obliquity promise to con- 
strain the planet's formation, tidal evolution, and tidal dis- 
sipation. What these data might reveal about a planet de- 
pends on whether the planet is unaffected by stellar tides, 
slightly aflfected by tides, or heavily influenced by tides. 

2.1. Tidally Unaffected 

The present-day rotation rate of a planet is the prod- 
uct of both the planet's formation and its subsequent evo- 
lution. Planets at sufficiently large distances from their 
parent stars are not significantly affected by stellar tides, 
thus these objects rotate with their primordial rates altered 
by planetary contraction and gravitational interactions be- 
tween the planet, its satellites, and other planets. To the 
degree that a planet's rotational angular momentum is pri- 
mordial, it may be diagnostic of the planet's formation. 
Planets formed in circular orbits from a protoplanetary 
disk inherit net prograde angular momentum from the ac- 
creting gas, resulting in rapid prograde rotation (Lissauer 
1995). Planets that form in eccentric orbits receive less 
prograde specific angular momentum than planets in cir- 
cular orbits, and as a result they rotate at rates varying 
from slow retrograde up to prograde rotations similar to 
those of circularly accreted planets (Lissauer et al. 1997). 
Thus, comparing the current-day rotation rates for planets 
in circular and eccentric orbits may reveal whether extra- 
solar planets formed in eccentric orbits or acquired their 
eccentricity later from dynamical interactions with other 
planets or a disk. 

The orientation of a planet's rotational axis relative to 
the vector perpendicular to the orbital plane, the planet's 
obliquity or axial tilt t, can also provide insight into the 
planet's formation mechanism. Jupiter's low obliquity 
(3.12°) has been suggested as evidence that its formation 




Fig. 3. — Anatomy of a transit, after BCGNB. 



was dominated by orderly gas flow rather than the stochas- 
tic impacts of accreting planetesimals (Lissauer 1993). 
Tidally unevolved extrasolar planets determined to have 
high obliquities could be inferred either to have formed 
differently than Jupiter or to have undergone large obliq- 
uity changes as has been suggested for Saturn {t = 26.73°; 
Ward & Hamilton 2002). 
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Fig. 4. — Oversimplified model of the effect of oblate- 
ness on a transit lightcurve. We plot the differences be- 
tween the lightcurve of an oblate planet and the lightcurve 
of a spherical planet with the same cross-sectional area, 
Ff=o,i{t) — Ff=o,o{t) while holding all other transit param- 
eters, R^, Rp, b, and ci the same and setting these values 
equal to those values measured by BCGNB for HD209458b. 
The top panel plots the lightcurve differential for impact 
parameters b ~ 0.0 (solid line), b = 0.3 (dashed line), and 
6 = 0.6 (dot-dashed line). An oblate planet for b < 0.7 
encounters first contact before, and second contact after, 
the equivalent spherical planet, resulting in the initial neg- 
ative turn for Ff=o.i{t)—Ff=o,o{t), and subsequent positive 
section of the curve for these impact parameters. In the 
bottom panel, we plot the differential lightcurve for b = 0.8 
(solid line), b = 0.9 (dashed line), and 6 = 1.0 (dot-dashed 
line). For b > 0.7, under otherwise identical conditions, 
an oblate planet first encounters the limb of the star af- 
ter the equivalent spherical planet, and last touches the 
limb later on ingress. For b = 0.9, because Rp > 0.1 i?* 
there is no second contact, i.e. the transit is partial, so the 
flux differential does not return to near zero, even at mid- 
transit. The middle plot shows the lightcurve differential 
at the changeover point between these two regimes, where 
6 = 0.7. 



2.2. Tidally Influenced 

Tidal interaction between planets and their parent stars 
slows the rotation of those planets with close-in orbits 
(Guillot et al. 1996). This tidal braking continues until 
the net tidal torque on the planet becomes zero. Whether 
a planet reaches this end state depends on its age, radius, 
semimajor axis, and the planet : star mass ratio. 

The rate of tidal braking also depends on the parameter 

Q, which represents the internal tidal dissipation within 
the planet. The value of Q is poorly constrained even for 
the planets in our own solar system (Goldreich & Soter 
1966). Nevertheless, measurements of extrasolar planet 
rotation rates could constrain Q for these planets based on 
the degree of tidal evolution that has taken place (Seager 
& Hui 2002). 

Tidal braking for objects with nonzero obliquity can, 
somewhat counterintuitively, act to increase an object's 
obliquity. Tidal torques reduce the component of a planet's 
angular momentum perpendicular to the orbital plane 
faster than they reduce the component of the planet's an- 
gular momentum in the orbital plane (Peale 1999). This 
occurs because at solstice the planet's induced tidal bulge 
is not carried away from the planet's orbital plane by plan- 
etary rotation. Therefore for large fractions of the year 
stellar tidal torques do not act to right the planet's spin 
axis, while torques that reduce the angular momentum per- 
pendicular to the orbital plane act year-round. As a result, 
planets that have undergone partial tidal evolution can ex- 
hibit temporarily increased obliquity as the planet's rota- 
tion rate decreases. Eventually, such a planet reaches max- 
imum obliquity and thereafter approaches synchronous ro- 
tation and zero obliquity simnultaneously. Planets under- 
going tidal evolution may be expected to have higher obliq- 
uities on average than planets retaining their primordial 
obliquity. 

2.3. Tidally Dominated 

The end state of tidal evolution for planets in circular 
orbits is a 1 : 1 spin-orbit synchronization between the 
planet's rotation and its orbital period, along with zero 
obliquity. However, most of the extrasolar planets discov- 
ered thus far are on eccentric orbits (Marcy et al. 2000) (we 
sometimes shorten 'planets on eccentric orbits' to 'eccen- 
tric planets' even though the eccentricity is not inherent to 
the planet). Thus these eccentric planets will never reach 
1 : 1 spin-orbit coupling as a result of tidal evolution be- 
cause the tidal torque (Eq. 2) on the planet from its star is 
much stronger (due to the dependence) near periapsis. 



The tidal torque between a planet and star is given by 
(Murray & Dermott 2000) 
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sgn(n - (j)) , 
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where k2 is the planet's Love number, R its radius, O its 
rotation rate (in radians per second), ^ the instantaneous 
orbital angular velocity (also in radians per second), and r 
is the instantaneous distance between the planet and star. 
The function sgn(a;) is equal to 1 if a; is positive and —1 if 
X is negative. The magnitude of the stellar perturbation is 
proportional to GM^, the product of the stellar mass (M*) 
squared and the gravitational constant (G). The planet is 
spiui down by tidal torques if its rotation is faster than its 
orbital motion {^l > (p), and it is spun up if its rotation is 
slower than the orbital motion < (p). 

If an eccentric planet were in a 1 : 1 spin-orbit state, 
it would be spun up by the star when its orbital angular 
velocity is greater than average near periapsis, and it would 
be spun down when the orbital angular velocity is low near 
apoapsis. The total positive tidal torque induced while the 
planet is in close will exceed the negative torque while it 
is far away, despite the shorter time spent near periapsis. 
Thus to be in rotational equilibrium with respect to stellar 
tides, an eccentric, fluid planet must rotate faster than its 
mean motion. 

The Earth's moon avoids this fate because it has a 
nonzero component of its quadrupole moment in its orbital 
plane. The torque that the Earth exerts on this permanent 
bulge exceeds the net tidal torque imparted on the moon 
due to its eccentric orbit, keeping the Moon in synchronous 
rotation (Murray & Dermott 2000). Fluid planets, how- 
ever, have no permanent quadrupole moment and thus 
have no restoring torque competing with the stellar tidal 
torques (Greenberg & Weidenschilling 1984). Tidal evolu- 
tion ceases for these bodies when the net torque per orbit 
is zero, which can only be achieved by supersynchronous 
rotation. 

The precise rotation rate necessary to balance the tidal 
torques over the eccentric orbit depends on how Q varies 
with the tidal forcing frequency (the difference between 
the rotation rate and instantaneous orbital angular veloc- 
ity, rip — 4>p). Conventionally, Q has been assumed to be 
either independent of the forcing frequency or inversely 
proportional to it (Goldreich & Peale 1968). The resulting 
equilibrium spin-orbit ratios as a function of eccentricity 
for each of these cases are plotted in Figure 1. 

Probably both of these assumptions for the behavior 
of Q are too simple. In particular, if the behavior of Q 
changes under a varying tidal forcing frequency (Op — cpp 
is a function of time), then the tidal equilibrium rotation 
rate would differ significantly from that plotted in Figure 
1. Measurement of rotational rates of eccentric extrasolar 
planets in tidal equilibrium could, in principle, either dif- 
ferentiate between these two models or suggest other fre- 
quency dependences, shedding light on the yet unknown 
mechanism for tidal dissipation within giant planets. 



3. ROTATION AND OBLATENESS 

Rotation affects the shape of a planet via two mecha- 
nisms: gravity must provide centripetal acceleration, thus 
the higher velocity at the equator causes the planet to 
bulge by transfer of mass from polar regions; and, secon- 
darily, the redistributed mass alters the planet's gravita- 
tional field and attracts even more mass toward the equa- 
torial plane. The ratio of the required centripetal accel- 
eration at the equator to the gravitational acceleration, q, 
represents the relative importance of the centripetal accel- 
eration term: 



Q = 



(3) 



where fl is the rotation rate in radians per second, Mp is 
the mass of the planet, and Rgq is the planet's equatorial 
radius (Murray & Dermott 2000). 

We use the Darwin-Radau relation to relate rotation 
and oblateness accounting for the gravitational pull of the 
shifted mass: (Murray & Dermott 2000): 
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where G is the planet's moment of inertia around the 
rotational axis and C is shorthand for GM~^R~^. The 
Darwin-Radau relation is exact for uniform density bod- 
ies (C = 0.4), but is only an approximation for gas giant 
planets (C ~ 0.25; Hubbard 1984). 

By combining Equation 3 and Equation 4, we arrive at 
a relation for rotation rate, O, as a function of oblateness, 
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For our solar system, the Darwin-Radau relation yields 
rotation periods accurate to within a few percent (Table 
1) using model-derived moments of inertia from Hubbard 
& Marley (1989). 

Extending the Darwin-Radau relation to extrasolar 
planets requires estimation of the appropriate moment 
of inertia coefficients, C, for those planets. For transit- 
ing planets whose masses and radii are known, we use a 
self-consistent hydrodynamic model of the planet and as- 
sumptions about its composition to estimate C following 
Fortney & Hubbard (2002). To first order, C is indepen- 
dent of oblateness due to similar symmetry around the 
rotational axis, so our hydrodynamic model does not need 
to explicitly incorporate oblateness. The Darwin-Radau 
relation and spherically symmetric hydrodynamic models 
provide sufficient precision for the current work; however, 
to estimate the oblateness as a function of rotation more 
robustly, a two-dimensional interior model involving both 
rotational and gravitational forces should be used (e.g., 
Hubbard & Marley 1989). 

Under the spherically symmetric assumption, we calcu- 
late the C of Juoiter to be 0.277 with no core and 0.253 



with a 20 core, and wc calculate the C of Saturn to be 
0.225 with a core (we are unable to calculate the interior 
structure of Saturn without a core due to deficiencies in our 
knowledge of the equation of state) . Using these moments 
of inertia instead of the measured ones listed in Table 1 
yields similar rotation errors of a few percent. We apply 
our model to generic 1.0 Rjup cxtrasolar giant planets of 
varying masses and architectures in Figure 2. 

For HD209458b, our models calculate C to be 0.218 with 
no core and 0.185 with a 20 core. To estimate the 
oblateness of HD209458b assuming synchronous rotation, 
we rearrange Equation 5 to solve for /, 
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and then use the synchronous rotation rate H, = 2.066 x 
10~^ radians per second to obtain / = 0.00285 with no 
core and / = 0.00256 with a 20 M0 core. These results 
imply an equator-to-pole radius difference of ~ 200 km, 
which, although small, is still comparable to the atmo- 
spheric scale height at 1 bar, ~ 700 km. 

Showman & Guillot (2002) suggest that zonal winds on 
HD209458b may operate at speeds up to ^ 2 km s^^ in the 
prograde direction, and Showman & Guillot (2002) go on 
to show that these winds might then spin up the planet's 
interior, possibly to commensurate speeds of several hun- 
dred m up to a few km (though the model of 
Showman & Guillot (2002) does not treat the outer lay- 
ers and interior sclf-consistcntly). This speed is a non- 
negligible fraction of the orbital velocity around the planet 
at the surface, 30 km s~-^, and is also comparable to the 
planet's rotational velocity at the equator, 2.0 km s^^. As 
such, if radiatively driven winds prove to be important on 
HD209458b they would affect the planet's oblateness. We 
use the rotation rate implied by the Showman & Guillot 
(2002) calculations to provide an upper limit for the ex- 
pected oblateness of HD209458b. If the entire planet were 
spinning at its synchronous rate plus 2 km s"'^ at the cloud 
tops, the rotational period would be halved to 1.8 days, 
with a corresponding oblateness of 0.0109 and 0.0098 for 
the no core and core models respectively. 

During revision of this paper, Konacki et al. (2003) an- 
nounced the discovery of a second transiting planet. This 
new planet, OGLE-TR-56b, has a radius of 1.3 Rjup, a 
mass of 0.9 Rjup, and an extremely short orbital period of 
1.2 days. Although these parameters are less constrained 
than those for HD209458b, wc proceed to calculate that 
the oblateness of this new object should be 0.017 with no 
core and 0.016 with a 20 Mgj core, for C of 0.228 and 0.204 
respectively. 

Current ground-based transit searches detect low- 
luminosity objects like brown dwarfs and low-mass stars 
in addition to planets. However, the high surface grav- 
ity for brown dwarfs and lower main sequence stars leads 
to low values of q and very small oblatenesses for those 
objects. For a 13Mjup brown dwarf with l.ORjup in an 
HD209458b-like 3.52 dav orbit, the expected oblateness is 



only 0.00007. Measuring oblateness will therefore only be 
practical for transiting planets and not for other transiting 
low-luminosity bodies. 

We also note that the measurement of oblateness along 
with an independent measurement of a planet's rotation 
rate ft would determine the planet's moment of inertia. 
This would provide a direct constraint on the planet's 

internal structure, possibly allowing inferences regarding 
the planet's bulk helium fraction and/or the presence of a 
rocky core. 

4. OBLATENESS AND TRANSIT LIGHTCURVES 
4.1. Transit Anatomy 

Brown et al. (2001) (hereafter BCGNB) investigated the 
detailed structure of a transit lightcurve while studying 
the Hubble Space Telescope lightcurve of the HD209458b 
transit. Figure 3 relates transit events to corresponding 
features in the lightcurve, modeled after BCGNB Figure 
4. The flux from the star begins to drop at the onset of 
transit, known as the first contact. As the planet's disk 
moves onto the star, the flux drops further, until at sec- 
ond contact the entire planet disk blocks starlight. Third 
contact is the equivalent of second contact during egress, 
and fourth contact marks the end of the transit. Due to 
stellar limb darkening the planet blocks a greater fraction 
of the star's light at mid-transit than at the second and 
third contacts, leading to curvature at the bottom of the 
transit lightcurve. 

For a given stellar mass, Af,, the total transit dura- 
tion, /, is a function of the transit chord length and the 
orbital velocity. We assume a circular orbit, which flxes 
the planet's orbital velocity. The chord length depends on 
the stellar radius, i?*, and the impact parameter, b. The 
impact parameter relates to i, the inclination of the orbital 
plane relative to the plane of the sky, as 



b = 



|acos(i) 



(7) 



where a is the semimajor axis of the planet's orbit. 



The duration of ingress and egress, w, is the time be- 
tween first and second contact, and is a function of Rp, 
i?*, and b (or i). Transits with 6 ~ have smaller w than 
transits with higher b. For transits with 6^1, called graz- 
ing transits, w is undefined because there is no second or 
third contact. 

The total transit depth, d, fixes the ratio Rp/R^, where 
Rp is the radius of the planet (except in the case of grazing 
transits). 

The magnitude of the curvature at the bottom of the 
transit lightcurve, r], determines the stellar limb darkening. 
We use limb darkening parameters ui and U2^ or ci and 
C2, which we define mathematically in Section 4.2. 



4.2. Methods 



We calculate theoretical transit lightcurves by compar- 
ing the amount of stellar flux blocked by the planet to the 
total stellar flux. The relative emission intensity across the 
disk of the star is greatest in the center and lowest along 
the edges as a result of limb darkening. Many parameteri- 
zations of limb darkening exist (see Claret 2000); however, 
we use the one proposed by BCGNB because it is the most 
appropriate for planetary transits. 

The emission intensity at a given point on the stel- 
lar disk, J, is parameterized as a function of /z = 
cos(sin~^(;9/i?*)), where p is the projected (apparent) dis- 
tance between the center of the star and the point in 
question. BCGNB defined a set of two limb darkening co- 
efiicients, which we call c\ and C2, that are are equivalent 
to 

(i-^)(2-^) (1-^)^ 
7(Yy = i-ci ^ + C2^— . (8) 

The advantage of this limb darkening function is that ci 
describes the magnitude of the darkening, while C2 is a cor- 
rection for curvature. This makes the BCGNB coefficients 
particularly useful for fitting transit observations because 
a good fit to data can be achieved by fitting only for the 
ci coefficient. 

Our algorithm for calculating the lightcurve takes ad- 
vantage of the symmetry inherent in the problem: that 
J(yu) depends only on p and not on the angle around the 
star's center, 6. We evaluate the apparent stellar flux at 
time r, F{t), relative to the out-of-transit flux F{), by sub- 
tracting the amount of stellar flux blocked by the planet 
from Fq: 
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(10) 
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where x(p, t) is the fraction of a ring of radius p and width 
dp covered by the planet at time r. In effect, we split 
the star up into infinitesimally small rings and add up the 
fluxes in Equation 9, then we determine how much of each 
of these rings is covered by the planet in Equation 10. 
We calculate the integrals numerically using Romberg's 
method (Press et al. 1992); x{r,t) is evaluated numeri- 
cally as well — there is no closed form general analytical 
solution to the intersection of an ellipse and a circle. 

This algorithm is more efficient than the raster method 
used by Hubbard et al. (2001) for planets treated as opaque 
disks because the use of symmetry and Romberg integra- 
tion minimize the number of computations of the stellar 
intensity, I{p). 
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Fig. 5. — Detectable difference between the lightcurve of 
an oblate (/ = 0.1) planet and the best-fit spherical model, 
fitting for R^,, Rp, 6, and ci. Higher b combined with 
larger values for iJ* and Rp simulate the lenghened ingress 
and egress of an oblate planet, diminishing the difference 
between oblate and spherical planets from Figure 4 for 
planets with b < 0.7 (upper panel: solid line, b = 0.0; 
dashed line, b = 0.3; dotted line, b = 0.5). For planets 
near b = 0.7, the length of ingress and egress cannot be 
simulated by higher b, and as a result the transit signal is 
highest for these planets (middle panel: solid line, b = 0.6; 
dashed line, b = 0.7; dotted line, b ~ 0.8). Above the 
critical value, b > 0.7, the oblate planet's signal can be 
simulated by lowering b for the spherical planet fit, reduc- 
ing the detectability of oblateness (lower panel: solid line, 
6 = 0.9; dashed line, 6 = 1.0). It is very difficult to deter- 
mine the oblateness for planets involved in grazing transits 
{b 1.0). The magnitude of the detectability difference is 
proportional to / to first order, hence to estimate the de- 
tectability of a planet with arbitrary oblateness, multiply 

the differences plotted here by g^;^ 
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4.3. Results 



To illustrate the effect oblateness has on a tran- 
sit lightcurve, we calculate the difference between the 
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Fig. 6. — Detectable difference between the lightcurve of 
an oblate (/ = 0.1) planet and its best-fit spherical model 
with 5 parameters: i?*, i?p, 6, ci, and C2. As in Figme 
5, differing values for i?*, Rp, and h for a spherical planet 
transit can allow it to mimic the transit of an oblate planet. 
With the addition of C2, the fit is much better for planets 
transiting at low impact parameter. (Upper panel: solid 
line, b = 0.0; dashed line, b = 0.3; dotted line, b = 0.5. 
Middle panel: solid line, b = 0.6; dashed line, b = 0.7; 
dotted line, b = 0.8. Lower panel: solid line, b — 0.9; 
dashed line, b = 1.0.) The differences for the 5 parameter, 

f 

like those for the 4 parameter model, vary as rpr — • 



lightcurve of an oblate, zero obliquity planet and that of a 
spherical planet with the same cross-sectional area. For fa- 
miliarity, we adopt the transit parameter values measured 
by BCGNB for the HD209458b transit: R.p = 1.347Rjup, 
= 1.146Ro, and Ci = 0.640. Plots of the oblate- 
spherical differential as a function of impact parameter 
are shown in Figure 4. 

For nearly central transits, an oblate planet encounters 
first contact before, and second contact after, the equiv- 
alent spherical planet. This situation causes the oblate 
planet's transit lightcurve initially to dip below that of the 
spherical planet. However, near the time when the planet 
center is covering the limb of the star, each planet blocks 
the same apparent stellar area and the stellar flux differ- 
ence is zero. As the two hypothetical transits approach sec- 
ond contact, this trend reverses and the spherical planet 
starts blocking more light than the oblate one until the 
oblate planet nearly catches up at second contact. Be- 
tween second and third contacts, the lightcurve differences 
are slightly nonzero because the two planets cover areas 
of the star with differing amounts of limb darkening. The 
lightcurves are symmetric, such that these effects repeat 
themselves in reverse upon egress. 

At high impact parameters (nearly grazing planet tran- 
sits) the opposite occurs. First contact for the oblate 
planet occurs after that for the spherical planet, because 
the point of first contact on the planet is closer to the pole 
than to the equator. In this scenario, the oblate planet's 
transit flux starts higher than, becomes equivalent to, and 
then drops below that of the spherical planet before return- 
ing to near zero for the times between second and third 
contact. In the case of a grazing transit, this sequence is 
truncated because there is no second or third contact. 

The boundary between these two regions occurs when 
the local oblate planet radius at the point of first contact is 
equal to i?ea, the radius of the equilvalent spherical planet. 
For planets that are small compared with the sizes of their 
stars {Rp <C R*) and that have low oblateness (/ < 0.1), 
this transition occurs when 9 = tt/A (from Figure 3) and 
b = \/2/2 = 0.707. The flux difference between transits of 
oblate planets and that for spherical planets, all else being 
equal, is at a minimum at this point and deviates from 
zero because the rate of change in stellar area covered is 
different for the two planets. First, second, third, and 
fourth contacts all occur at nearly the same time for each 
planet. However, if all else is not equal, as is usually the 
case since the stellar parameters are poorly constrained, 
then this result is misleading and does not represent the 
detectability of planetary oblatness. 



5. TRANSIT LIGHTCURVES AND EXOPLAN- 
ETS 



To test whether these flux differences are detectable, 
we use a Levenberg-Marquardt curve-fitting algorithm 
adapted from Press et al. (1992) to fit model transits 
to both the HST HD209458b lightcurve and hypothetical 
model-generated lightcurves by minimizing y^. As a test. 
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Fig. 7. — Best-fit parameters resulting from fitting a simu- 
lated transit liglitcurve for a planet with oblateness / using 
a spherical planet transit model. Due to the degeneracy in 
measuring transit parameters i?*, Rp, and b, and oblate- 
ness /, the best spherical fit to the oblate data have larger 
radii and higher impact parameters than the actual val- 
ues (see Section 5.1). Objects with negative values of / 
are prolate, a physically unreasonable proposition that we 
include here for completeness. 
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Fig. 8. — Here we plot the change of fitted spherical model 
parameters with oblateness as a function of impact pa- 
rameter. The derivitive of impact parameter with respect 
to oblateness continues off the top of this graph to 3.0 at 
b = 0.0. An oblate planet with zero obliquity transiting 
near the center of the star induces much higher deviations 
of the fitted parameters from the actual parameters than 
would a planet transiting near the critical impact param- 
eter, b = 0.707. 



i?* = 1.142Rq, Rp = 1.343Rjup, 
Cl = 0.647, and C2 = —0.065 with a reduced 
consistent with the values obtained by BCGNB. 



86.72° (6 = 0.504), 
" of 1.05, 



In order for planetary oblateness to have a noticeable 
effect on a transit lightcurve, it must be distinguishable 
from the lightcurve of a spherical planet. For a spherical 
planet transit model, the combination of transit parame- 
ters that correspond to the lightcurve that best simulates 
the data from an actual oblate planet transit become the 
measured values, and these measured quantities may not 
be similar to the actual values. Therefore to consider the 
detectability of planetary oblateness, we compare oblate 
planet transit lightcurves to those of the best-fit spherical 
planet lightcurve instead of to the lightcurve of a planet 
that differs from the actual values only in the oblateness 
parameter (as we did in Section 4.3). 

5.1. Zero Obliquity 

We compare a model transit of an oblate planet (/ = 
0.1) with zero obliquity, as is the case for a tidally evolved 
planet, to the transit of the best-fit spherical planet in Fig- 
ure 5 and Figure 6. The oblate planet transit signature 
is muted in each when compared to Figure 4 due to a de- 
generacy between the fitted parameters i?*, Rp, b, and the 
oblateness /. This degeneracy is introduced by the uncon- 
strained nature of the problem: in essence, we are trying 
to solve for 5 free parameters, i?*, Rp, 6, ci, and /, given 
just 4 constraints, d, w, I, and rj assuming (as BCGNB 
did) knowledge of the stellar mass Af* . Without assuming 
a value for , absolute timescales for the problem vanish, 
yielding a similar conundrum of solving for Rp/R^,, b, ci, 
and / given only d, rj, and w/l. The previous two sentences 
are intended only to be simplifications, as at higher photo- 
metric precision more information about the conditions of 
the transit is available. Hereafter we assume knowledge of 
M* , though this analysis could also have been done without 
this assumption, or with an assumed relation between M, 
and i?* as proposed by Cody & Sasselov (2002). Changes 
in i?*, Rp, and b mimic the signal of an oblate planet by 
altering the ingress and egress times while maintaining the 
total transit duration by keeping the chord length con- 
stant. 

For planets transiting at low impact parameter {b < 
0.7), an oblate planet's longer ingress and egress (higher w 
from Figure 3) are fit better using a spherical model with 
a higher impact parameter than actual, thus lengthening 
the time between first and second contact. Since for a 
given star transits at higher impact parameter have shorter 
overall duration, the best-fit spherical model has a larger 
i?* than the model used to generate the data to maintain 
the chord length, and thus a larger Rp to maintain the 
overall transit depth. 

Similarly, for simulated lightcurve data from a transit- 
ing oblate planet at high impact parameter {b > 0.7), the 
best-fit spherical model has a lower impact parameter than 
the simulated planet to increase the duration of ingress and 
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Fig. 9. Dctcctability of oblatcncss and obliquity relative 
to the best-fit spherieal model for planets with nonzero 
obliquity. The top panel shows the detectability differenee 
for b = 0.35 and t = 0° (dotted line), t = 15° (dot-dashed 
line), t = 30° (dashed hue), and t = 45° (solid line). The 
middle panel represents the same varying obliquities cal- 
culated for b = 0.7, and the bottom panel shows the dif- 
ferences for b = 0.83. The shape of the difference is quali- 
tatively similar for each case. The difference is maximized 
for t — 45° and b = 0.7, and falls off as the parameters 
near t = 0°, t = 90°, b = 0.0, and b = 1.0. Due to the 
symmetry of the problem, obliquities the ones shown plus 
90° are the time inverse of the differences shown. 



egress. For these planets, the fitted spherical parameters 
have smaller i?, and Rp than the simulated planet to main- 
tain the character of the rest of the lightcurve. 

Oblate planets that have impact parameters near the 
critical value, b ~ 0.7, cannot be as easily fit using a spher- 
ical model because changes in the impact parameter of 
the fit cannot increase the duration of ingress and egress. 
For these planets, the difference between the oblate planet 
transit lightcurve and the best-fit spherical planet tran- 
sit lightcurve is maximized, providing the largest possible 
photometric signal with which to measure oblateness. 

At present, it is necessary to fit for i?* and stellar limb- 
darkening parameters because of our limited knowledge of 
these values for the host stars of transiting planets. If 
i?* were known to sufficient accuracy (less than 1%), then 
with knowledge of the degeneracy between i?*, Rp, 
and b would be broken, allowing measurement of planetary 
oblateness without fitting. However, without assumptions 
about the stellar mass, knowledge of i?, would only help 
to constrain M,. Cody & Sasselov (2002) show that, for 
the star HD209458, current evolutionary models combined 
with transit data serve to measure the stellar radius to a 
precision of only 10%. 

In Figure 5 we plot the difference between the transit 
lightcurve of a hypothetical planet with the characteristics 
of HD209458b but an oblateness / = 0.1, and that planet's 
best-fit spherical model fitting for R^,, Rp, b, and ci. For 
low values of the impact parameter, b < 0.5, the best-fit 
spherical lightcurve emulates the oblate planet's ingress 
and egress while leaving a subtly different transit bottom 
due to the planet traversing a differently limb-darkened 
chord across the star. The magnitude of the difference 
is approximately a factor of 10 smaller than the non-fit 
di£Fc;rc;nce from Figure 4. Near the critical impact param- 
eter, b = 0.7, the transit lightcurve bottoms are very sim- 
ilar since the best-fit b is very similar to the actual oblate 
planet's impact parameter; however, the ingress and egress 
differ in flux by a few parts in 10^ for / = 0.1. For grazing 
transits, 6 ~ 1.0, the best-fit spherical model's lightcurve is 
indistinguishable from the oblate planet transit lightcurve 
even at the 10~^ relative accuracy level. 

Figure 6 shows the same differences as in Figure 5, but 
with a second stellar limb darkening parameter, BCGNB's 
C2, also left as a free parameter in the fit. Including C2 in 
the fit allows the fitting algorithm to match the transit bot- 
tom (the time between second and third contacts) better. 
This effect leads to excellent fits of oblate planet transits 
using spherical planet models and reduces the detectable 
difference to less than one part in 10^ for 0.0 < 6 < 0.5 and 
b > 0.9 using HD209458 stellar and planetary radii and an 
oblateness of / = 0.1. For 6 = 0.0 and 6 = 1.0 the spheri- 
cal model emulates an oblate transit particularly well, with 
the differences between the two being only a few millionths 
of the stellar flux. 

In both the 4- and 5-parameter zero-obliquity models, it 
is easiest to measure the effects of oblateness on the transit 
light curve for transits near the critical impact parameter. 



For HD209458 values with / = 0.1, the oblatcncss signal 
then approaches 3 x 10~^ for tens of minutes during ingress 
and egress and peaks near b = 0.8 instead of the critical 
value of 6 = 0.7 due to the finite radius of HD209458b 
relative to its parent star. 

Based on observations of the Sun, Borucki et al. (1997) 
expect the intrinsic stellar photometric variability on tran- 
sit timescales to be ~ 10~^. Jenkins (2002) used 5 years of 
3 minute resolution SOHO spacecraft data to deduce the 
noise power spectrum of the Sun. These noise effects are 
close enough in magnitude to the transit signal so as to 
affect the dctectability of a transiting oblate planet. 

Future high-precision space-based photometry missions 
such as MOST and Kepler may be able to detect the ef- 
fect for highly oblate transiting extrasolar planets, but the 
S/N ratio could be so low as to make unambiguous mea- 
surements of oblateness diffucult to obtain. 

If an observer were to fit photometric timeseries data 
from a transit of an oblate planet without fitting explicitly 
for the oblateness / (as, for instance, if the precision is 
insufficient to measure /), the planet's oblateness will be 
manifest as an astrophysical source of systematic error in 
the determination of the other transit parameters. Figure 
7 shows the effect that oblateness has on the stellar radius, 
planetary radius, impact parameter, and limb darkening 
coefficient for HD209458b. This variation is a strong func- 
tion of the planet's impact parameter. In Figure 8, we show 
the how this systematic variation changes as a function of 
the initial impact parameter for the HD209458b system. 
As a severe but still physical example, an IID209458b-like 
planet with oblateness / = 0.1 that transits at 6 = 0.0 
would be measured to have radii 2% above actual and an 
impact parameter nearly 0.3 above the real impact param- 
eter. 

5.1.1. HD209458b 

To illustrate the robustness of the degeneracy between 
i?*, Rp, b, and / discussed in Section 5.1, we fit the HST 
lightcurve from BCGNB using a planet with a fixed, large 
oblateness of 0.3(!). The best-fit parameters were i?* = 
1.08 Rq, Rp = 1.26 Rjup, b = 0.39, and ci = 0.633 with 
a reduced = 1-06 — indistinguishable in significance 
from the spherical planet model! Although unlikely, a high 
actual oblateness for HD209458b could alter the measured 
value for the planet's radius into better agreement with 
theoretical models. 

The expected dctectability of oblateness for HD209458b 
is extremely low. During ingress and egress the transit 
lightcurve for IID209458b should differ from that of the 
best-fit spherical model by only one part in 10^ for / = 0.01 
and at the level of 3 x 10^^ if / = 0.003. For comparison, 
the BCGNB HST photometric precision is 1 x 10"^. 

Although the systematic error in the determination of 
transit parameters can be important for highly oblate plan- 
ets, it is not at all significant for HD209458b. Assuming 
HD209458b has an oblateness of / = 0.003 as calculated in 



Section 3, the fitted radii are only ^ 0.05% above the ac- 
tual radii, and the fitted impact parameter is 0.0006 above 
what the actual impact parameter should be. 

5.1.2. OGLE-TR-56b 

Measuring the oblateness of the new transiting planet 
OGLE-TR-56b (Konacki et al. 2003), / = 0.016 (see Sec- 
tion 3), would require photometric precision down to at 
least 4 X 10~^. Since the impact parameter for this object 
is as yet unconstrained, the above precision corresponds 
io b — 0.7. For other transit geometries, higher precision 
photometry would be necessary to measure the oblateness 
of this object. 

5.2. Nonzero Obliquity 

We plot the dctectability of oblateness and obliquity 
for planets with nonzero obliquity in Figure 9. Here we 
define the projected obliquity (which we refer to as just 
obliquity hereafter), or axial tilt t, as the angle between the 
orbit angular momentum vector and the rotational angular 
momentum vector projected into the sky plane, measured 
clockwise from the angular momentum vector (see Figure 
3). The major effect of nonzero obliquity is to introduce 
an asymmetry into the transit lightcurve (Hui & Seager 
2002). 

The plots in Figure 9 are difference plots for planets with 
different obliquities and impact parameters, yet all have 
the same shape qualitatively. In trying to fit the asym- 
metric ingress and egress lightcurves, the best-fit spherical 
planet splits the difference between them. For planets with 
< t < 7r/2, this process yields a difference curve that 
initially increases as the spherical planet covers the star 
at a faster rate than does the oblate planet. Because of 
the asymmetric nature of the problem, however, the egress 
of the oblate planet takes longer than the spherical one, 
causing an initial upturn near third contact. The transit 
lightcurve difference for planets with > t > — 7r/2 is equal 
to the one for < t < 7r/2, but reversed in time. 

This general shape is the same for the asymmetric com- 
ponent of each transit lightcurve, and it is superimposed 
on top of the symmetric component studied in Section 5.1. 
The asymmetric component is maximized near the critical 
impact parameter (6 = 0.7), because the planet crosses the 
stellar limb with its projected velocity vector at an angle of 
7r/4 with respect to the limb. For transits across the mid- 
dle of the star, 6 = 0, the asymmetric component of the 
lightcurve vanishes as a new symmetry around the planet's 
path is introduced, and the local angle between the veloc- 
ity vector and the limb is 7r/2. Similarly the asymmetric 
planet signal formally goes to zero for grazing transits at 
b = 1.0, but in practice the asymmetric component is still 
high for Jupiter-sized planets. 

Planets with obliquities of zero {t = 0) are symmetric 
and have no asymmetric lightcurve component (see Section 
5.1). Likewise, planets with t = tt/2 are also symmetric. 
The asymmetric component is maximized for planets with 



t = 7r/4. 

Transiting planets with nonzero obliquity c;an break the 
degeneracy between transit parameters discussed in Sec- 
tion 5.1. However, when the obliquity is nonzero a de- 
generacy between projected obliquity and oblateness is in- 
troduced: a measured asymmetric lightcurve component 
of a given magnitude could be due to a planet with low 
oblateness but near the maximum detectability obliquity 
oit = 7r/4, or it could be the result of a more highly oblate 
planet with a lower obliquity. In this case only a lower limit 
to the oblateness can be determined based on the oblate- 
ness for an assumed obliquity of 7r/4. This degeneracy can 
be broken with measurements precise enough to determine 
the symmcrtrie lightcurve component. In addition, transit 
photometry is only able to measure the projected oblate- 
ness and obliquity of such objects due to the unknown 
component of obliquity along the line-of-sight (Hui & Sea- 
ger 2002), therefore the true oblateness is never smaller 
than the measured, projected oblateness. 

The asymmetric transit signal of an oblate planet with 
nonzero obliquity could also be muddled by the presence 
of other bodies orbiting the planet. Orbiting satellites or 
rings could both introduce asymmetries into the transit 
lightcurve that may not be easily distinguishable from the 
asymmetry resulting from the oblate planet. Satellites 
around tidally evolved planets are not stable (Barnes & 
O'Brien 2002), and rings around these objects may prove 
to be difficult to sustain as well. However, objects that are 
not tidally evolved, those farther away from their parent 
stars, may potentially retain such adornments, and their 
effects on a transit lightcurve could be difficult to differen- 
tiate from oblateness. 

Orbital eccentricity can also cause a transit lightcurve to 
be asymmetric. Although the eccentricity of HD209458b 
is zero to within measurement uncertainty based on radial 
velocity measurements, future planets discovered solely by 
their transits may not have constrained orbital parameters. 
For these objects with unknown eccentricity, if the eccen- 
tricity is very high then under some conditions the velocity 
change between ingress and egress may be high enough to 
emulate the oblateness asymmetry discussed in this sec- 
tion. However, we do not explicitly treat that situation in 
this paper. 

6. CONCLUSIONS 

Examining a transiting planet's precise lightcurve can 
allow the measurement of the planet's oblateness and, 
therefore, rotation rate, beginning the process of charac- 
terizing extrasolar planets. To a reasonably close approxi- 
mation (a few percent), the rotation rate of an extrasolar 
giant planets is related to the planet's oblateness by the 
Darwin-Radau relation. Measurements of a planet's rota- 
tion rate could constrain the tidal dissipation factor Q for 
those planets, as well as possibly shed light on the tidal 
dissipation mechanism within giant planets based on the 
spin : orbit ratios of tidally evolved eccentric planets. 



The detectable effect of oblateness on the lightcurve of 
a planet with zero obliquity is at best on the level of a few 
parts in 10^. This effect may be discernable from space 
with ultra high precision photometry, but could prove to 
be indistinguishable from stellar noise or other confound- 
ing effects. Accurate independent measurements of stellar 
radius can break the degeneracy between the stellar radius, 
planetary radius, and impact parameter, allowing for much 
easier measurement of oblateness. Without such measure- 
ments, the primary effect of oblateness and zero obliquity 
when studying transit lightcurves will be to provide an as- 
trophysical source of systematic error in the measurement 
of transit parameters. 

Planets with nonzero obliquity have higher detectabil- 
ities (up to ~ 10"'*) than planets with no obliquity, but 
yield nonunique determinations of obliquity and oblate- 
ness. In order to obtain unique obliquities and oblatenesses 
for these objects, precision comparable to that needed for 
the zero obliquity case is needed. 

The detectability of oblateness for transiting planets is 
maximized for impact parameters near the critical impact 
parameter of = 0.7. Many transit searches are currently 
underway, yielding the potential for the discovery of many 
transiting planets in coming years. Given the opportu- 
nity to attempt to measure oblateness, our analysis sug- 
gests that the optimal observational target selection strat- 
egy would be to observe planets around bright stars that 
transit near an impact parameter of 0.7. 

The authors wish to acknowledge Bill Hubbard and Fred 
Ciesla for useful conversations; Wayne Barnes, Robert 
H. Brown, Christian Schaller, and Paul Withers for 
manuscript suggestions; and Robert H. Brown for valu- 
ble advice. 
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Table 1 

Darwin-Radau in the Solar System 



Planet 


f 


C 


Calculated Period 


Actual Period 


Error 


Jupiter 


0.06487 


0.264011 


10.1609 hr 


9.92425 hr 


2.38% 


Saturn 


0.09796 


0.22037^ 


10.8817 hr 


10.6562 hr 


2.12% 


Uranus 


0.02293 


0.2268 1 


16.6459 hr 


17.24 hr 


3.45% 


Neptune 


0.01708 


0.23 * 


16.8656 hr 


16.11 hr 


4.69% 


Earth 


0.00335 


0.3308 


23.8808 hr 


23.9342 hr 


0.223% 



^Hubbard & Marley (1989) 
* Assumed 

Note. — Comparison of calculated rotation rates from the Darwin-Radau ap- 
proximation (Equation 4) to actual rotation rates for selected planets in our solar 
system. 



